Attractive electron-electron interaction induced by geometric phase in a Bloch band 
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We investigate electron pairing in the presence of the Berry curvature field that ubiquitously exists 
in ferromagnetic metals with spin-orbit coupling. We show that a sufficiently strong Berry curvature 
field on the Fermi surface can transform a repulsive interaction between electrons into an attractive 
one in the p-wave channel. We also reveal a topological possibility for turning an attractive s-wave 
interaction into one in the p-wave channel, even if the Berry curvature field only exists inside the 
Fermi surface (circle). We speculate that these novel mechanism might be relevant to the recently 
discovered ferromagnetic superconductors such as UGe2 and URhGe. 



Attractive interaction between electrons (or neutral 
fermions) is responsible for superconductivity (or super- 
fluidity). In condensed matter systems, attractive inter- 
action is usually induced by the boson-exchange mech- 
anism pj]. Indeed, in the celebrated Bardeen-Coopcr- 
Schreiffer (BCS) theory, electrons develop attractive in- 
teraction by exchanging phonons Subsequent stud- 
ies show that other collective excitations such as charge 
density waves Q and spin fluctuations [4] can also induce 
attractive interaction. The boson-exchange mechanism, 
together with the concept of Cooper-pair, is considered 
to be the cornerstone of the modern theory of super- 
conductivity, and its validity can even be extended to 
the unconventional superconductors such as the high- Tc 
cuprates 

In this Letter, we show a new possibility for the occur- 
rence of attractive electron-electron (e-e) interaction in 
ferromagnetic metals with spin-orbit coupling. We take 
the ferromagnetic state as given, and focus on the effect 
of the Berry curvature field which exists ubiquitously in 
such materials [sl, 0] • Our question is relevant because 
superconductivity has been found within ferromagnetic 
phase, such as in UGe2 and URhGe [§]. The Berry 
curvature effect on electron motion is analogous to a mag- 
netic field in the reciprocal space 13, 11 1 , and has been 
invoked to successfully explain the anomalous Hall effect 
in ferromagnets [12, 13, 14]. On the other hand, unlike 
a magnetic field in real space, monopole sources for the 
Berry curvature field can occur in the reciprocal space at 
band degeneracy points. In the vicinity of the monopoles, 
the Berry curvature becomes very strong. 

Our theory is formulated within an effective one band 
model, where ferromagnetism and spin-orbit coupling 
has already been taken into account, such as one cal- 
culated self-consistently from a spin-density functional 
theory. Figure [ija) shows the origin of the Berry curva- 
ture field: an electron evolving adiabatically in the re- 
ciprocal space will accumulate a geometric (Berry) phase 
4>B = j^{uk\idkUk) ■ dfe associating with the adiabatic 
change of the quasi-momentum k [ist , in analogy to the 
Aharanov-Bohm phase acquired by electron moving in 
the real space in the presence of a magnetic field. It 




Figure 1: (a) Electron moving adiabatically in the reciprocal 
space acquires a Berry phase; (b) Berry curvature field in 
the vicinity of a band degeneracy split by magnetization and 
spin-orbit coupling can be modeled as a "magnetic field" in 
the reciprocal space generated by a "monopole" out of the 2D 
Brillouin manifold. 



suggests a fictitious "magnetic field" in the reciprocal 
space with the "vector potential" A.{k) — {uk\idkUk) 
and the corresponding "physical field" (Berry curvature 
field) rj(fc) = Vfc X A{k), where Uk is the periodic part of 
the Bloch wave function for the electron band concerned. 

The central result of this work is that attractive inter- 
actions in the p-wave channel may be produced with the 
help of the Berry curvature field. We show that the pres- 
ence of a sufficiently strong Berry curvature field on the 
Fermi surface can transform a repulsive e-e interaction 
into an attractive one in the p-wave channel. There is 
also a topological effect analogous to the Aharanov-Bohm 
phase. This is for a situation where the Berry curvature 
field vanishes or is negligible on the Fermi surface but 
not so inside of it. A Berry phase around the Fermi sur- 
face can still result from the flux within. We show that 
an originally attractive interaction in the s-wave channel 
can be turned into one in the p-wave channel. 

To be specific, we investigate the following effective 
one-band many-body Hamiltonian: 



H = J2e{k,)+J2vin-f,), 



(1) 



where e(fc) is the quasi-particle dispersion operator, and 
V{ri — rj) is the two-body e-e interaction, k (r) is 
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the quasi-momentum (position) operator, and the in- 
dexes i, j denote the particle number. By assuming the 
usual canonical commutation relations for k and r, the 
Hamiltonian had served as the basis for many successful 
condensed matter theories, for instance, the conventional 
BCS theory. 

Our study is motivated by the understanding that 
in the presence of the Berry curvature field, the differ- 
ent components of the position operator r do not com- 
mute [iE[illI3: 



(2) 



where = x,y,z denote the different components of 
the position operator. It is then interesting to see how 
the change of the electron dynamics dictated by the non- 
commutative position operator affects the electron corre- 
lations in our system 

To proceed, we introduce the canonical coordinates R 
which satisfy the usual commutation relations [i?^, R^] = 
and [Rfj,, k^] — i(5^,y, which can be realized if we define: 



R = r-Aik), 



(3) 



where A.{k) = {uk\idkUk) is the "vector potential" 
corresponding to the Berry curvature field n. We 
can then second-quantize the many-body Hamiltonian 
Eq. 11]). We first re-express the interaction poten- 
tial in terms of the canonical coordinates R. This 
can be done by making use of the Fourier expansion 
V{f, - Tj) = {27r)-'^ J dqv{q)e"i-^^e-"i-^^ and the re- 
lation e"i-^ = e'xiqMe^Q-Re-^xiq.k)^ where x is defined 
by the equation q ■ Vfcx(q, k) = q ■ ^(fe). The interac- 
tion potential can then be expressed in the plane wave 
basis ipkiR) = 1/ VV exp{ik ■ R), from which its second- 
quantized form can be easily read out, yielding finally: 



H = J2<l^)cick 



1 

2V 



E 

kk'K 



u{k,k';K)c 



■-+k' 



(4) 



where V is the total volume of the system, (c) is the 
quasi-particle creation (annihilation) operator, and 



u{k,k';K) =v{k' -k)e 



i4>B{f+k',^+k)+i4,Bif-k',^-k) 

(5) 

i.e., the interaction is modified by a geometric phase de- 
fined as (f>B{k, k') ~ A.{k) -dk with the integral along 
the straight line connecting k and k' in the reciprocal 
space. (j}B is exactly the Berry phase acquired by an 
electron scattered from k to k' . In Eq. (|4]), we omit the 
spin index, and focus on ferromagnetic systems in which 
the spin degrees of freedom are fully quenched due to the 
strong magnetization and spin-orbit coupling. 

In the following, we investigate how the geometric 
phase modifies the e-e interaction. We notice that the 



distribution of the Berry curvature field in a Brillouin 
zone is governed by the fc-points of band degenera- 
cies, which in mathematics are equivalent to "magnetic 
monopoles" in the reciprocal space 13, 15, T^. We thus 
focus on one of such "monopoles" and see how the geo- 
metric phase in its vicinity modifies e-e interaction. To 
identify the essential physics without being obscured by 
complexities in mathematics, we limit our study on a two 
dimensional (2D) ferromagnetic system. In such a sys- 
tem, there is usually no band degeneracy However, 
one can usually find band near-degeneracies at high sym- 
metry points of the Brillouin zone which are only split by 
the presence of magnetization and spin-orbit coupling. 
In the vicinity of these points, the Berry curvature field 
can be modeled as a "magnetic field" in the reciprocal 
space generated by a "monopole" out of the 2D Bril- 
louin manifold, as shown in Fig. [^b). We thus have: 
n{k) = (OM/2)KB/(fc2 + 4)3/2fc, and 



Aik) 




Kb 



(6) 



where we assume that the "monopole" is located at 
(0,0, Kb). Qm — ±1 is the charge of the "monopole". 
Kb measures how close the 2D system is to the band 
degeneracy. It is related to the magnitude of the band 
gap induced by the magnetization and spin-orbit cou- 
pling: the larger band gap, the larger kb (for instance, 
see Ref. [H). 

Attractive e-e interaction induced by the Berry curva- 
ture field: First, we demonstrate that the strong Berry 
curvature field in the vicinity of a reciprocal space "mag- 
netic monopole" could transform a repulsive e-e inter- 
action to an attractive one in p-wave channel. The 
e-e interaction in a typical metal can be modeled as 
Vir) — Vo exp[— ktf(v''^^ + a^^a)]/^(r/a)2 + 1, which 
is screened at large distances for r ^ 1/ktf and satu- 
rates at small distances for r a, where ktf is the 
Thomas- Fermi screening wave vector Its Fourier 

transformation reads. 



v{q) = vq- 



exp 



■irpp — KTF 



vW^tf) 



1 



(7) 



with vq — 2T:Voa/ Ktf- Equation ([7]) should be consid- 
ered as the renormalized interaction between "dressed" 
electrons resulting from a complete treatment of a bare 
many-body Hamiltonian [21]. 

For the isotropic model considered here, the effective 
interaction between a pair of electrons with the opposite 
momentum {K — 0) can be classified by their relative an- 
gular momentum = mh. For channel m, the effective 
potential is 2l|: 



1 r'^'^ 

^{k, k')^— dev{k' - fc)e2'*B(fe''fe)e»™«, (8) 
Jo 



Figure 2: Effective e-e interaction Um- (a) Typical behav- 
ior of Um(k,k) for different m. Parameters: ^sa = 0.2 and 
ktfIi^b ~ 2.6; (b) fc-position of the Um=i{k, k) minimum as 
a function of ktf/ks for different values of kbo, from 1 to 10; 
(c) The minimum (most attractive) value of the effective po- 
tential for p-wave channel (m = 1) as a function of ktf/i^b 
for different values of kbci- The different curves are offset 
vertically for clarity. The dotted line shows the boundary de- 
termined from Eq. for the onset of attractive interaction. 
Qm = 1. 



where 9 is the angle from k to k' and we have made use 
of the relation <j>B{—k', — fc) = 0_b(A;', k). 

Figure mja) shows the the effective interaction for dif- 
ferent channels for a given set of parameters. Attrac- 
tive interaction (i.e., Um{k, k) < 0) is evident for channel 
m — I (p-wave) and rn = 3 (/-wave). Unlike the con- 
ventional attractive interaction due to boson exchange 
which is always present in a thin shell near the Fermi 
surface [ij, the effective interaction due to the geometric 
phase is attractive only in the vicinity of the "monopole" 
(i.e., k = 0). Figure. [2l^b-c) show the dependence of the 
effective interaction (fc-position where the effective inter- 
action is the most attractive, and its magnitude, respec- 
tively) on the parameters {kb, i^tf, a) for the p-wave 
channel (m = 1). We note that the attractive interaction 
only occurs in certain regime of the parameter space. 

The condition for the onset of attractive effective in- 
teraction can be determined by examining the limit of 
fc, fc' ^ kb,ktf, where v{q) « exp [-q^/2Kl] with 
K„ = ktf/V^ + KTFd and the Berry curvature field can 
be considered as a constant flzik) ~ ^Iq = Qa//2k^, 
with the corresponding geometric phase: 



(PB{k,k') = -nokk' sine. (9) 



Figure 3: Phase diagram in ep (Fermi energy) -T (Temper- 
ature) plane, es = e(Ks). Filled dots indicate the region 
occupied by the superconducting phase. Solid line shows 
the superconducting gap Af(0) at the Fermi surface and 
at the zero temperature, scaled by a factor 1/1.76. The 
good correspondence between the phase boundary and the 
solid line suggests the usual BCS relation AF{0)/kTc ~ 
1.76. The dashed line shows an empirical fitting Af(0) ~ 
aes exp[—l / pF\um{kF , fe^)!], where pp is the density of states 
at Fermi surface and a 5. The reasonably good fitting sug- 
gests the correlation between the magnitude of the supercon- 
ducting gap and the strength of the attractive interaction at 
the Fermi surface. Parameters: KBa — 0.2, htf/hb = 2.6, 
and pfVq = 2. The electron dispersion is assumed to be the 
simple parabolic form. Qm = 1. 



Then it follows that 



Um (fcj fc ) 





1-0O 


m/2 


■ fc2-hfc'2- 






l + fZ-n 


exp 


[ 2^1 \ 


X 



[ {~irJ.r (^sgn(0o)y0^^) , \M > 1 ' 

where 0n = i^o^^- Jm {Im) is the (modified) Bessel 
function of the first kind. The attractive interaction (i.e., 
Um{k,k) < 0) arises in the channels with odd positive 
(negative) m for JIq > (flo < 0) if 

\no\nl>l. (11) 

The boundary determined from this criteria is shown as 
the dotted line in Fig. Wl,c) , which coincides well with the 
boundary directly determined from the numerical result. 

We apply the BCS gap equation in Ref. [2l| to in- 
vestigate the superconducting phase induced by the at- 
tractive interaction. The result is summarized in Fig. [S] 
The superconducting state has p-wave symmetry with 
A(A;) = A{k){kx ± iky) (for Qm = ±1)- The magnitude 
of the superconducting gap strongly depends on the po- 
sition of the Fermi surface, which is a direct result of the 
strong fc-dependence of the effective potential. 
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The superconductivity we predict is closely associ- 
ated with ferromagnetism, which breaks the time-reversal 
symmetry, and together with spin-orbit coupling, gives 
rise to the Berry curvature field in the vicinity of the 
high symmetry fe-points. In this picture, the super- 
conducting phase naturally coexists with ferromagnetism 
and disappears when the ferromagnetism is suppressed. 
This behavior makes it a plausible alternative theory 
for the recently discovered ferromagnetic superconduc- 
tors UGe2 H and URhGe % In the traditional pic- 
ture, enhanced spin fluctuations near a quantum criti- 
cal point are responsible for the pairing of electrons. It 
predicts the superconducting phase on both sides of the 
ferromagnetic-paramagnetic transition point [j^ . which 
contradicts with the experimental finding that the su- 
perconducting phase only exists in the ferromagnetic 
side. On the other hand, we note that the conditions 
for the onset of superconductivity with the mechanism 
(i.e., Eq. (fTTjl and the Fermi surface must reside in the 
vicinity of the "monopole" ) is rather stringent . Further 
investigations is required for establishing the definite con- 
nection between the theory and real systems. 

Unconventional pairing symmetry of the topological 
origin: Second, we consider the case that the Fermi- 
surface (-circle) is far from the "monopole", i.e., kp ^ 
Kb- While there is no strong presence of the Berry cur- 
vature field on Fermi-circle in this case, the "monopole" 
still presents a reciprocal space "magnetic flux" thread- 
ing through the Fermi-disc with a total flux — t^Qm- 
The corresponding geometric phase in the vicinity of the 
Fermi-circle reads: 



6B(fc,fc')«%^ 



(12) 



Using Eq. 
is: 



the effective e-e interaction in channel m 



i{k,k') « Vm+QM{k,k') 



(13) 



where Vm is the Fourier component of the bare e-e inter- 
action at channel m. For an originally attractive inter- 
action in s-wave channel (i.e., Vm=o < 0), the effective 
interaction Um is attractive in channel m = —Qm — 
giving rise the p-wave pairing symmetry. The uncon- 
ventional pairing symmetry is of the purely topological 
origin: were the "magnetic flux" not present, the bare 
interaction would favor the s-wave pairing symmetry. 

The mechanism presents a new possibility distinctly 
different from the conventional scenario in which the e-e 
interaction induced by the Boson-exchange mechanism 
has the particular structures (e.g., hot spots induced 
by Fermi-surface nesting) that favors the unconventional 
pairing symmetry. In contrast, the unconventional pair- 
ing symmetry here is a result of the intrinsic structure of 
the Bloch band, originated from the "monopole" (band 
degeneracy) buried deep inside the Fermi-sea. 



So far, our discussion is built upon the Berry cur- 
vature field associating with the adiabatic evolution of 
quasi-electrons in reciprocal space. The approach makes 
possible to construct a minimum theory for the effects 
of the geometric phase in a Bloch band to electron cor- 
relations. On the other hand, it relies on the assump- 
tion that the interaction potential varies slowly over the 
atomic length scale. When the potential is not slowly 
varying, we have to work with the matrix elements be- 
tween the Bloch states, il'ki'r) = e^'^"^\uk{r)). Ignoring 
inter-band terms and Umklapp contributions, as usually 
justified, we obtain an effective one-band Hamiltonian 
just like Eq. g]) but with (for K = Q) 

u{k, k';K = 0) = v{k' - k){uk'\uk){u-k'\u-k)- (14) 

In this case, the extra phases are now the Pancharatnam 
geometric phases 23|. For smooth potentials, only for- 
ward scattering is important for which fc' — fc is small, 
where the Pancharatnam phases reduce to the Berry 
phases we discussed before [2J]. 
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